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Abstract 

We address the persistence of Holder continuity for weak solutions of the linear drift-diffusion equation 
with nonlocal pressure 

ut + b • Vu - At* = Vp, V • u = 

on [0, oo) x R™, with n > 2. The drift velocity b is assumed to be at the critical regularity level, with respect 
to the natural scaling of the equations. The proof draws on Campanato's characterization of Holder spaces, and 
uses a maximum-principle-type argument by which we control the growth in time of certain local averages of 
u. We provide an estimate that does not depend on any local smallness condition on the vector field 6, but only 
on scale invariant quantities. 

1 Introduction 

A classical problem in partial differential equations is to address the regularity of solutions to parabolic problems 
involving advection by a vector field b and diffusion 

u t + b- Vu - Ait = 0. (1.1) 

If the vector field b is sufficiently regular, the solution u is expected to be regular as well. Naturally, this is 
expressed as a result of the type: if b is bounded with respect to some norm, then u is smooth in some sense. The 
appropriate norms for such statement are the ones that are either critical or subcritical with respect to the inherent 
scaling of the equation. More precisely, if u is a solution of dl.lt . then for any r > the function u r (x, t) = 
ru(rx, r 2 t) solves an equation of the same form but with drift velocity given by b r (x, t) = rb(rx, r 2 t). This 
change of variables acts as a zoom in that focuses on the local behavior of the solution u. An assumption on 
b is critical with respect to the scaling of the equation if the norm of b r coincides with the norm of 6, for any 
r > 0. The assumption would be subcritical if b r has smaller norm than b for all small enough values of r, and 
supercritical otherwise. 

As a rule of thumb, with current methods it seems impossible to obtain a regularity result for (ll.lt with a 
supercritical assumption on b, since the transport part of the equation would be stronger than diffusion in the 
small scales. With subcritical assumptions on b, it is generally possible to treat equation (ll.lt as a perturbation 
of the heat equation, and strong regularity results in this direction are available. The Kato class condition for 
b is probably the largest class that falls into this category. For results in the subcritical case see for example 

mama. 

Obtaining regularity estimates for (ll.lt depending only on scale invariant norms of b requires the use of 
non-perturbative techniques, since the drift term does not become negligible at any scale. To the best of our 
knowledge, the only results available are variations of the De Giorgi-Nash-Moser Harnack inequality lfl8l . which 
states that weak solutions to ( ll.lt are C a for positive time, for some small a > 0. Results in this direction 
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include a variety of critical assumptions on the vector field b. For b G I%L% with n/q + 2/p ~ 1, we refer 
to lfl6l Chapter 3] or |[T9l . For divergence-free drift b G LfBMO^. 1 , the Holder regularity of weak solutions 
was proved only recently in (9) and 11221 . For b in a space-time Morrey space, this result was obtained recently 
in |[T9l . See ifTTl |2"T1 |26l for other conditions on b yielding Holder regularity, such as the form boundedness 
condition. 

The equation (II. lb is essentially a scalar equation, since even if u is a vector field, each component would 
satisfy the same equation. In contrast, the De Giorgi-Nash-Moser theory is hard to apply to actual systems. In 
this article we consider a Stokes system with drift, i.e. we add a pressure term as is common for the equations 
of fluid dynamics, and we look for a solution u which is divergence free. Given a divergence free vector field 
b : [0,T] x R" — > M. n , we consider the following evolution equation 

u t + b ■ Vu - Ait = Vp (1.2) 

for a solution u : [0, T] x M" — > M" which satisfies 

V-u = 0. (1.3) 

The term pressure gradient may be computed from ( |1.2) - (|1.3) by the formula 

Vp = V(-A)" 1 div(6- Vm). (1.4) 

In this paper we prove that if a scale invariant norm of b is bounded, then the C a norm of u at time t is 
bounded in terms of its C a norm at time zero. Our result is a propagation of regularity instead of a regularization 
result, in the sense that we require the initial data uq to be Holder continuous, and a G (0, 1) is arbitrary. 

The assumption on the divergence-free drift velocity b is that it is an L p integrable function in time, with 
values in the L 1 -based Morrey-Campanato space M 13 , where j3 G [—1, 1], and p = 2/(1 + fi). We recall 
cf. 11241 Definition 1.7.2] the definition of the i 1 -based Morrey-Campanato spaces M s . For any s > —n, let 
N = max{-l, [s]}. We say / G M s if / e L\ and 



sup sup r s inf — / \f(z) - P(z)\ dz < +oo, (1.5) 
xei»o<Ki -PeP" \B r (x)\ J Br {x) 

where V N is the set of all polynomials of degree less than or equal to N, with the convention that V^ 1 = {0}. 
In this article, we will only consider f3 in the range [—1, 1]. The condition b(-, t) € A/* 3 has a different character 
depending of the value of the exponent /?, cf. Il24l Theorem 5.3.1]: if (3 = 1, M@ coincides with the space of 
Lipschitz functions (here we subtract elements of V instead of V x , the later yielding the Zygmund class) ; if 
f3 G (0, 1), it is exactly the Holder class C"; if j3 = 0, it corresponds to the class of functions of bounded mean 
oscillation BMO; while if /3 £ [—1, 0) it is the usual Morrey-Campanato space. In all these cases, the estimate 
in our main theorem depends only on the semi-norm [&(•, t)] M fi associated to the space. In this paper we consider 
divergence-free drifts b such that 

[b{;t)] MP = sup swpr-^M(x,t,r) < g(t) (1.6) 

kGR 3 >">0 

for some g G L p ([0,T]), where we define 



M(x,t t r) = \ [ \b(z,t)-b(x,r,t)\ dz= [ \b(x + ry,t) -5(z,r,*)| dy (1.7) 

I B r (x) JBi(O) 



jtfi 



and b(x, r, t) is chosen to equal zero if /3 € [—1, 0), the average of b over B r {x) if (3 = 0, respectively b(x, t) if 
f3 G (0, 1], which is equivalent to ( 11.51 ) (except for (3 = 1). We give further details on the precise assumptions 
on b in Section [3] below. Our main theorem in the case /3 G ( — 1, 1] is given in Theorem 11.11 below, while the 
endpoint case /? = — 1 is addressed in Theorem |4.1| (see also Remark [OV 
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Theorem 1.1. Assume b : [0, T] x W n — > M" /s a divergence-free vector field such that b £ £ p ([0, T]; M' 3 ) w'f/z 
/3 G (— 1, 1] and p = 2/(1 + /3). Asrame a/so f/taf wo € C a for some some a £ (0, 1). 77?e« f/tere ex/sfs a weaA; 
solution u : [0, T] x R" — > K" of the system 

ut+b-Vu- Au = Vp (1.8) 
V-w = (1.9) 
u(x,0)=u (x) (1.10) 

5Mc/i f/taf u(x, t) is C a in xfor all positive time t £ (0, T]. Moreover, we have the estimate 

M;t)]o < C[uo] c -, 
for some positive universal constant C = C(T, a, (3, [^Il^m 13 )■ 

Remark 1.2. To the best of our knowledge, the result of Theorem 11. 11 is new even if the pressure term was 
removed, and the scalar equation (11. It was considered instead. Indeed, the De Giorgi-Nash-Moser iteration 
scheme provides a C a estimate for the solution for some small a only, whereas our result provides a C a estimate 
for any a £ (0, 1). 

Remark 1.3. The assumption b £ L^M^ implies a local smallness condition in the sense that |H|_Lp([t-T, 
becomes arbitrarily small as r — > 0, due to uniform integrability, but without any rate. This is not true for the 
endpoint case p = oo. However, this is not the reason why we require p < oo in Theorem ll.il and in fact this 
local smallness plays no role in our proof. Indeed, the constant C in the estimate of Theorem 1 1 . 1 1 depends only 
on the scaling-invariant norm of b, and not on any other feature of the vector field b, such as the modulus of 
continuity of the map Q H> II ^|| z, p j\^' 3 (Q) ■ argument that relies on the local smallness of b would make the 
constants in the estimates depend on the rate at which the local norm of b decays, and would hence be implicitly 
a subcritical result. 

In the particular case when b — u, Theorem 11.11 becomes a no-blowup condition for solutions to Navier- 
Stokes equation. It says that if the norm of u remains bounded in L\MP, then u does not blow up on R 3 x [0, T]. 
This is a scale invariant condition that is slightly more general (but has the same scaling) than the classical 
Ladyzhenskaya-Foias-Prodi-Serrin condition u G L v t L q x with 2/p + n/q < 1 (cf. m[l5]|20]|23)). N °te that the 
endpoint case (q,p) = (3, oo), when n = 3, was only treated recently in ifTUl (see also ifTTl and Theorem 14.11 
below for a related statement). The result of Theorem 1 1.1 1 for the full range of p 6 [1, oo) may be new for the 
Navier-Stokes equations as well, though it is comparable to other available regularity criteria in terms of scaling 
critical norms of u (cf. [5] |6] [14) and references therein). 

One important difficulty for proving Theorem II. H is to deal with the non-local pressure term on the right 
hand side of (ll.2l i. There are very few results of this kind available for equations with pressure terms. In ||261 the 
same equation ( 11.2b is considered and a Lipschitz estimate is shown under a sub-critical assumption on b (which 
includes b £ L^L% with 2/p + n/q < 1 — e, for any e > 0). 

The idea of the proof is to write the Holder regularity condition of u(-, t) in integral form using a classical 
theorem of Campanato [2). Then we claim that these local integral estimates have a certain growth in time (in 
terms of integral estimates on b). In order to prove that these estimates hold for all time we argue by contradiction 
and look for the first point in which they would be invalidated. At that time we apply the equation and obtain 
a contradiction in a way that resembles maximum-principle-type arguments (see also lfl2l [131 for the SQG 
equations). The integral representation of the Holder modulus of continuity allows us to take advantage of the 
divergence-free condition and the integral bound on b. The divergence free condition on u is used in the estimate 
for the gradient of pressure term. The general method of the proof introduced here seems to be new, as it may 
be applied to systems with pressure gradients, and we believe it may be applicable to other evolution equations 
in the future. 

In section [4] we analyze the endpoint case /3 = — 1. The method of this article is applicable in this case, but 
we need to impose an extra smallness condition on the vector field b (cf. (14. U - fi^l i below). 
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We believe that the most important contribution of this article is the introduction of a new method to prove 
Holder estimates for evolution equations. Nevertheless, we show an example of how Theorem 11.11 can help 
prove that a nonlinear equation is well posed. Let us consider the following modified energy critical Navier- 
Stokes equation in 3D (see Q for a similar modified critically dissipative SQG equation) 

d t u+ (i-A)- 1/4 u-V^u- Au = Vp (1.11) 
divu = (1.12) 

that is, b = (- A)" 1/4 u in dOt-dOl). It follows directly from TheoremfTTlthat the system ( 11.1 lb — ( II. 12b is well- 
posed in the classical sense. Indeed, the global existence of weak solutions u £ Lf 3 L 2 r\L 2 H 1 , is straightforward, 
as is the local existence of strong solutions. In particular, for any T > 0, we have that {—A)~ x l A u is a priori 
bounded in L 2 (0, T; iJ 3 / 2 (R 3 )), and by the Calderon-Zygmund theorem we have 



<oo (1.13) 

L 2 (0,T;BMO) 



for any T > 0. Therefore, by applying Theorem 1 1.1 1 with (3 = 0, we obtain that u(-, t) £ C a on (0, T], given 
that u(-, 0) = uo £ C a , where a £ (0, 1) is arbitrary. From this estimate, it is easy to obtain higher regularity 
of u by a standard bootstrap argument. 

Note that for the above system one could also use classical energy estimates at the level of vorticity, combined 
with Sobolev interpolation, to obtain the global well-posedness of the problem. On the other hand our method 
allows some extra flexibility in the relationship between b and u. As explained above, when b = {— A)^ 1 ' 4 u 
the system is well posed. Following essentially the same idea we can obtain using Theorem 1 1.1 1 that the system 
is well posed for any of the following choices 

• b = a(x)( — A) _1 / 4 u + Vg for any bounded function a in R 3 and \7q is the gradient of a scalar function 
that makes b divergence free. In this case we apply the a priori estimate u £ L 4 L 3 that is obtained by 
interpolation from the energy inequality, and gives b £ L 4 L 6 . 

• b = J k(x, y)u(y) dy where k(x, — ) £ £, 2n /i n + 2 ) f or anv x g jj™ an( j k(—, y) is divergence free for any 
y £ W l . 

• Ab = div(u <g> u) + Vq. 

We plan to explore other applications of this method in the future. 



2 Preliminaries 

In this section we state a few introductory remarks about the weak and classical solutions to ( 1 1 .8b — ( fTT9b . and 
recall a classical characterization of Holder spaces in terms of local averages. Throughout the rest of the paper 
we will write WW to denote L p t = LP(0, T; L q ), and similarly LPM? will be used instead of L v t M%. 

We first prove Theorem 11.11 assuming that the solution is classical (i.e. C 2 in space and C 1 in time). The 
important feature is that the a priori estimate (11.6) depends only on the assumptions of Theorem [TTT] and not on 
any further smoothness assumptions on b or u. Then we approximate any weak solution with classical solutions 
by using a mollification of b, and pass to the limit to obtain the result of Theorem ll.ll in full generality. 

Definition 2.1 (Weak Solutions). Ifb £ £^([0, T] x R") is divergence-free, a function u £ L°°([0, T] x R") is 
a weak solution of d 1 -2 b ). if it is weakly divergence-free, and for all smooth, divergence-free, compactly supported 
test functions (p we have: 

/ tp(x,T)u(x,T) dx + / u (— (ft + bV(p — Aip) dx dt = / Lp(x, 0)u(x, 0) dx. 

JR" i[0,T]xR™ JR" 
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The following proposition is standard. 

Proposition 2.2. Let b e and u e be a sequence of smooth divergence-free vector fields. Assume that u e is a weak 
solution of ( 11.2b with drift velocity b e . Assume also that b e — > b in L} oc Lui oc Then, up to a subsequence, u £ 
converges weakly to a weak solution of M.2\ . 

Using proposition ^. 21 we immediately observe the following. 

Proposition 2.3. It is enough to prove Theorem U.l \ assuming that b is smooth and u is a classical solution. 

Proof of Proposition ^. 3\ The assumption b E LPM^ implies in particular that b E L\ oc l} uloc . Using a mol- 
lification argument, we consider a sequence of smooth vector fields b e converging strongly to b in Lj^L^ . 
Moreover, we choose b e such that ||fr £ Hlpmi 3 i s bounded uniformly with respect to e (for example mollifying 
with a smooth function with fixed L 1 norm). For each of these vector fields, we solve the equation (11.21 ), for 
instance using the mild formulation and Picard iteration, to obtain a smooth solution u e . If the result of Theo- 
rem ll.ll is known for classical solutions, we would have that u e satisfies the estimate ( 11.6b uniformly in e. Note 
that in particular we also obtain u e £ L°° . By Proposition ^. 21 up to a subsequence, u e converges weakly to a 
weak solution u of jl.2\ , and therefore this solution u satisfies dl.61 l as well. □ 

In order to prove the main theorem, we use a local integral characterization of Holder spaces. For this 
purpose, let (p be a nonnegative, radially symmetric, smooth function supported in Bi(0). Unless otherwise 
specified, the center of the unit ball B\ in M" shall be 0. Let us also assume that J ip(y) dy = 1. The following 
theorem (or a small variation of it) is proved in |]2]- 

Theorem 2.4 (Campanato's characterization of Holder spaces). Let f : R" — s- K m be an L 2 function such 
that for all r > and x £ R n , there exists a constant f such that 

f \f(x + ry)~f\ 2 p(y)dy<A 2 r 2a (2.1) 

for some positive constant A, and a € (0, 1). Then the function f has a Holder continuous representative such 
that 

\f(x)-f(y)\<BA\x-y\ a 
where the constant B depends on dimension and a only. 

The most natural choice of the constant / in the above theorem, for which the converse also holds, is to 
choose the average of / in the ball 

/ = / f(x + ry)<p(y) dy. 

This is optimal in the sense that it minimizes the left hand side in (12. It (see also dl.5t ). 

The theorem of Campanato is interesting because it provides a non-obvious equivalence between a Holder 
modulus of continuity, which is a priori a pointwise property, and averages of differences of the function, which 
is an integral property. This relation will allow us to exploit the divergence free nature of the vector fields b and 
u when estimating the evolution of a Holder modulus of continuity. 

3 Evolution of a modulus of continuity 

We will prove that the solutions of ( 11-8 b — ( fTTTOb do not lose regularity by showing that they always satisfy a 
time dependent Holder modulus of continuity. This modulus of continuity will evolve and deteriorate with time, 
but it will stay bounded. In order to take advantage appropriately of the divergence-free character of the vector 
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field u, we use the integral characterization of the modulus of continuity. Let ip be a radially symmetric weight 
supported in B\ with mass one as in section[2] We denote the weighted mean of u on B r (x) by 

u(x,t,r) = u(x + ry,t)ip{y) dy. (3.1) 

The integral version of the modulus of the continuity of u is then 

I(x,t,r) = \u(x + ry,t) - u(x,t,r)\ 2 ip(y) dy. (3.2) 

Due to Theorem 12. 41 if we knew that 

I{x,t,r) <f(t) 2 r 2a , (3.3) 

for some function f(t) > 0, and all r > 0, then [u(-, t)]o < C fit) for some universal constant C. Our goal 
is to prove that estimate ( 13.3b holds for all t > 0, if it holds at t = 0, for some function f(t) to be chosen 
appropriately. 

As discussed in the introduction, our assumptions on b will be in terms of quantities similar to 1, which are 
distinguished by the parameter /3 £ [—1, 1] as follows. 

(i) The Morrey-Campanaoto case. For f3 £ [—1,0), let 

M(x,t,r)= [ \b(x + ry,t)\dy = \ [ \b(z,t)\ dz. (3.4) 

JBx T JBr-{x) 

We assume that there exists a positive function g £ £ / 2 ^ 1+ ^ sucn that 

sup svpr- p M(x,r,t) < g(t) ||6(-,t)||^ < g(t) (3.5) 

i6t" r>0 

for all t > 0, where || • \\ m s denotes the usual Morrey norm (cf. Il24l ). 

(ii) The BMO case. For (3 = 0, we let 

M{x,t,r)= [ \b(x + ry,t)-b{x,r,t)\dy = \[ \b(z,t) -b(x,r,t)\ dz, (3.6) 
Jb x r n J Bt[x) 

b(x,r,t) = \ I b(z,t) dz (3.7) 

r JB r (x) 

is the usual mean of b on B r (x). We assume that there exists a positive function g £ L 2 such that 

sup supM{0,x,r,t) < g(t) ||6(., t)||sA/o < (3-8) 

for all < > 0, where || ■ || bmo denotes the norm on the space of functions with bounded mean oscillation. 

(iii) The Holder and Lipschitz cases. For j3 £ (0, 1], we consider 

M(x,t,r)= [ \b{x + ry,t)-b(x,t)\dy = \ \b{z,t) - b{x,t)\ dz. (3.9) 

JB 1 r JB r (x) 

We assume that there exists a positive function g £ L 2 such that 

sup sup r-PM(x,t,r) < g(t) & [b(;t)] c o < g(t) (3.10) 

igt" r>0 

for all t > 0, where [-] c p denotes the Holder semi-norm. Note that 2/(1 + f3) £ [1, 2) when (3 £ (0, 1]. 



where 
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We shall prove that if ( 13.5b . ( 13.8b . or respectively (13.101 holds, then we have I(x, t, r) < f(t) 2 r 2a for all 
t > 0. The proof is in the flavor of a maximum principle. We show that if the inequality is satisfied at t = 0, it 
will be satisfied for all positive t. For the critically dissipative SQG equation, which as opposed to ( 11.8 b — ( fTT9b is 
a scalar equation, and has an L°° maximum principle, Kiselev, Nazarov, and Volberg fT3l , use an argument in 
the same spirit, but where the breakdown is considered for the pointwise modulus of continuity. 

In order to prove Theorem 11.11 assume there is a first time t and some value of x where the strict modulus is 
invalidated, i.e. 

l(x,t,r) = f(t) 2 r 2a . (3.11) 

By Proposition 12. 31 we can assume that u is a smooth function vanishing at infinity. Therefore the equality in 
( 13. 1 II ) of the modulus must be achieved at some r > and x <G K™. 

If we fix t and r, the function / achieves its maximum at x, and we obtain 

= V X 1 = f (u(x + ry,t)-u(x,t,r))-(W x u{x + ry,t)-W x u(x,t,r))(p(y)dy. (3.12) 

Due to the definition of u (13.11 ), and the fact that V x u(x, t, r) does not depend on y, we also have 

0=/ (u(x + ry, t) - u(x, t, r)) • V x u(x + ry, t)ip(y) dy. (3.13) 

JBi 

Since I < f(t) 2 r 2a for all times prior to t, and since [u — u)ip has zero mean, we thus conclude 

2f'(t)f(t)r 2a < d t I = [ {u(x + ry,t)-u{x,t,r))-{d t u(x + ry,t)-d t u(x,t,r))<p(y)dy (3.14) 

J Bi 

(u(x + ry, t) - u(x, t, r)) ■ d f u(x + ry, t)ip(y) dy. (3.15) 

The key to prove Theorem ll.ll is to find an appropriate upper bound for the right hand side of (13.151 ) in terms 
of f(t) and A4(r, t). Inserting the equation (11 .2b in the right hand side of ( 13.151 ), we obtain 

2f'(t)f(t)r 2a < f (u(x + ry,t)-u(x,t,r)) 

JBt 

■ ( - b(x + ry, t) ■ V x u{x + ry, t) + A x u(x + ry, t) - Vp{x + ry, t)) tp(y) dy 

= A + V + V. (3.16) 

The following three lemmas give bounds to the three terms on the right side of ( 13.16b . The advection term A is 
the simplest one to estimate. Observe that 

V x u(x + ry) — -V y u(x + ry). (3.17) 

r 

This identity, together with the assumption div 6 = 0, allows us to integrate by parts the gradient into the weight 
(p and obtain a precise estimate for A. 

The dissipative term turns out to be negative, but we must analyze it with care in order to obtain a precise 
lower bound on its absolute value. In fact, note that if u is linear in B r (x) then T> = 0. We will obtain an 
estimate of T> that measures how much u is forced to separate from a linear function, just from the values of I 
and d r X at the point where the equality ( 13.1 lb holds. 

Lastly, we obtain an upper bound for the pressure term V, comparable to the advection term A. This is to be 
expected since Vp is obtained from b ■ Wu though an operator of order zero. However, the pressure estimate is 
more involved since the formula for the pressure is non-local and in order to obtain the desired estimate we need 
to take advantage of some cancellations that occur after integration by parts of Riesz kernels using that both b 
and u are divergence free. 

We now carry out the estimates for the three terms on the right of ( 13.161 ) in the three lemmas below. 
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Lemma 3.1 (The advection term). Let u and b be as in the statement ofTheorem \l.l\ Then we have 

A = [ (u(x + ry, t) - u(x, t, r)) ■ (-b(x + ry, t) ■ V x u{x + ry))<p{y) dy < Cr 2 ^ 1 f(t) 2 M(x, t, r) 

(3.18) 

for all /3 G [—1, 1], where C is a positive constant depending only on a, ip, and n. 

Proof of Lemma \3J\ Using the identities (13.13b , (13.17b , and integrating by parts, we obtain from ( 13.18b that 

A=- I \u(x + ry,t)-u(x,t,r)\ 2 (b(x + ry,t)-b(x,r,t))-Vip(y)dy (3.19) 



where b = b(x, r, t) is a constant with respect to y, to be chosen suitably in the three cases for f3 € [—1,1], as 
discussed above. From identity ( 13.191 1, the Holder inequality, and Theorem |2.4| we directly obtain 

A^Cr^^fitfMix^^), (3.20) 

for all j3 € [—1, 1], where C is a positive constant depending on a, n, and tp through sup Bi |Vy|. □ 

The second term corresponds to the viscosity and it is strictly negative, as we will show below. 

Lemma 3.2 (The dissipative term). Let u and b be as in the statement of Theorem U .1\ Then we have 

V= { A x u(x + ry,t)-(u(x + ry,t)~u(x,t,r))(p(y)dy<-cf(t) 2 r 2a - 2 (3.21) 



for all r > 0, where c is a sufficiently small positive constant, depending only on n, a, and tp. 



Remark 3.3. Note that the constant c in Lemma l3T2| goes to zero as a —> 1. This is the reason why Theorem 
1 1.1 1 works for a < 1 only. 



In order to prove Lemma 1372] we need the following technical result, which relates the quantities V and I. 
Lemma 3.4. For a fixed x, we have 

c mr)-^r)f (3 . 22) 

for some sufficiently small positive constant c. 

For simplicity we omited the x dependence in d3.221 i. Note that if / is a linear function, then T> = and 
I = Cr 2 for some constant C. Moreover, in this case rl' = 21. We see that if (rl' — 21) is non zero, then the 
function u cannot be linear. The inequality ( 13.221 i gives a precise quantitative version of this fact. 

We note that the identity 



\f(y)-c\ 2 p(y)dy= / / (f(y) - f {z)) 2 p{y)p{z) dy dz (3.23) 

JBi JB X 

holds for any constant c, in particular for c = /, giving equivalent definitions for I and T> in terms of double- 
integrals. The proof of identity (13.23b is straightforward, and hence omitted. 

With formula ( 13.23b in mind, we prove the following lemma, which is exactly the case r = 1 for d3.221 i. The 
proof of Lemma l3~4l follows for all other values of r from Lemma l3~31 bv scaling. 

Lemma 3.5. Let f : B± — > R™ be any H 1 function. There is a constant C depending only on n and ip such that 
(f(y) - f(*))MvM*) dydz- [[ (y Vf(y) - z ■ V/(z)) ■ (f(y) - f(z))cp(y)cp(z) dy dz 



B 1 xB 1 JJB 1 xB 1 

1/2 , » ~ \ 1/2 



<C( // (Vf( y )-Vf(z)) 2 p(y)<p(z)dydz\ // {f {y) - f {z)) 2 V {y) V {z) dy dz 

3 1 xB 1 / \JJb 1 xB 1 J 
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Proof of Lemma \J75\ We start by writing /(y) — f(z) as an integral of V/ along the segment between z and y. 
Thus, we have 



// (/(v)-/(*))V(wM*) 

J JB 1 xB 1 



dy dz 



(V/(sy + (1 - s)z) ■ (y - z)) ■ (f(y) - f(z))<p(y)<p(z) ds dy d 



z. 



'B 1 xB 1 JO 

Substituting in the first term of the left hand side, we obtain 



(f(y) - f(z))MvM*) d?/ dz 



2 / / (y V/(y) - « ■ V/(z)) • (/(y) - f{z))<p(yMz) dy dz 



((V/(sy + (1 - s)z) - V/(y)) • y - (V/(sy + (1 - - V/(y)) ■ z) 
x (/(y) - /MMyVO) ds dy dz 

< I i/a f i/a j 
where I is as defined in ( 13.2b . and £ is defined as 



£ 



ixBi JO 



((V/(sy + (1 - s)z) - V/(y)) ■ y - (V/(sy + (1 - s)z) - V/(y)) • z)' ds dy dz. 



The quantity £ is a bounded quadratic functional with respect to the vector field V/ that vanishes whenever V/ 
is constant. Therefore, there is a constant C such that 



£<C (V/(y) - V/(z))V(yMz) dy dz = CV, 

J J B t xBi 

which concludes the proof. □ 

Proof of Lemma \J?2\ We recall that at the first point of equality in (13.1 It . the function I achieves its maximum 
as a function of x. Then we have A X I < 0, and hence 

> / (A x u(x + ry, t) - A x u(x, t, r)) • (u(x + ry, t) - u(x, t, r))kp{y) dy 

JBi 

Therefore, since (u — u)tp has zero mean and A x u does not depend on y, we obtain the inequality 

£><-/ \Vu(x + ry,t)-V x u(x,t,r)\ 2 ip(y)dy. (3.24) 

J Bi 

According to (13.24b , T> is negative, but we must estimate how negative it is. It is convenient to rewrite the 
formula for I{x, t, r) using (13.23b as 

X(x,t,r)= I \u(x + ry,t)-u(x,t,r)\ 2 <p(y) dy 

\u(x + ry, t) — u(x + rz, t)\ 2 ip(y)(p(z) dy dz. (3.25) 

Bi xBi 
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At the point where f(t) 2 r 2a = X(x, t, r) for the first time, we have S7 X X = and A X X < 0, so that 

> A x X(x,t,r) 



= 2 / / \V x u(x + ry,t)- V x u(x + rz,t)\ 2 

J J B 1 xB 1 

+ (A x u(x + ry, t) — A x u(x + rz, t)) ■ (u(x + ry, t) — u(x + rz, t))ip(y)ip(z) dy dz 
= 2V + 2 / / \V x u(x + ry, t) - V x u(x + rz, t)\ 2 y(y)(p(z) dy dz, 

J JB 1 xB 1 

where T> is given by ( 13.241 i. We thus have proven the inequality 

T><- \V x u(x + ry,t) - V x u(x + rz,t)\ 2 ip(y)tp(z) dy dz 

J J Bi xfli 

= — j // \V y u(x + ry,t)-V z u{x + rz,t)\ 2 Lp(y)^(z)dydz. (3.26) 

The right hand side is clearly negative unless u is an affine function. We now need to estimate how negative 
it is, in terms of I and d r I, Note that at the first point of equality 1 = f(t) 2 r 2a , we must also have d r X = 
2af(t) 2 r 2a ~ 1 . We compute <9 r Ias 



d r I = 2 j (u(x + ry, t) - u(x, t, r)) y ■ V x u(x + ry, t)ip(y) dy 

(u(x + ry, t) — u(x + rz, t)) ■ (y ■ \7 x u(x + ry, t) — z ■ \7 x u(x + rz, t))tp(y)ip(z) dy dz 



'B 1 xB 1 

(u(x + ry, t) - u(x + rz, t)) ■ (y ■ V y u(x + ry,t) - z ■ V z u(x + rz, t))ip(y)ip(z) dy dz. 

(3.27) 

From the expression (13.25) and (I3.27K we can apply Lemma 1331 to obtain 

(22 - rd r l) < Cl 1/2 ([ [ \X7 y u(x + ry,t) -V z u(x + rz,t)\ 2 (p(y)tp(z) dy dz) 

Recalling the inequality ( 13.261 ), we obtain 

V < -4 (2J ~: W = -c/(t) 2 (l - a) V- 2 (3.28) 
r l X 

for some positive constant c, depending only on a, n, and (p. □ 

Lastly, we bound the pressure term V arising on the right side of (13.161 . 
Lemma 3.6 (The pressure term). Let u and b be as in the statement ofTheorem \l.l\ Then we have 

[u(x + ry, t) - u(x, r, t)) ■ Vp(x + ry, t)<p{y) dy < C/(t) 2 5 (i)r 2 «+^ 1 (3.29) 



for all r > and f3 G [—1,1], where C is a positive constant, depending on a, /3, n, and (p. 



Proof of Lemma UM Since all estimates in this section hold for a fixed time t > 0, we omit the time dependence 
of all functions. Recall that the function b = b(x, r) (which is constant respect to y) is chosen to be b = b(x) in 
the Holder case, the average of b over B r (x) in the BMO case, or b = in the Morrey-Campanato case. 
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In order to estimate the third term in (13.16) . let us analyze the identity ( 11.4b . Since u is divergence free, we 
have 

V x p(x + ry) = V a; (-A) 2 : 1 div x (b(x + ry) ■ V x u(x + ry)) 

= -Vyi-Ay)- 1 dW y ((b(x + ry) - b(x , r)) ■ V„(u(3; + ry) - u(x , r))) . (3.30) 

We have the pressure term equal to 



(u(x + ry) - u(x, r)) ■ V x p(x + ry)tp(y) dy 
= i j ' (u(x + ry) ~ u(x, r)) ■ Vy(-Ay) -1 div y ((b(x + ry) - b(x, r)) ■ V y (u(x + ry) - u(x, r))j <p(y) dy. 
We integrate the V y by parts and use that u is divergence free to obtain that V equals 

- J (-Ay) -1 div y (jb(x + ry) - b(x, rj) ■ V y {u{x + ry) - u(x, r))j (u(x + ry) - u(x, rj) ■ Vf(y) dy. 
Using that b and u are divergence free, we re-write the above identity as 

V = - J (-Ay^didj ((b(x + ry) - b(x,r))i(u(x + ry) - u(x, r))^ (u(x + ry) - u(x)) ■ V<p(y) dy 
= j J RiRj ((b(x + ry) - b(x,r))i(u(x + ry) - u(x , r)) j^j (u(x + ry) - u(x,r)) ■ V<p y (y) dy 
= — j (b(x + ry) — b(x 7 r))i(u(x + ry) ~ u(x,r))jRiRj(^u(x + ry) — u(x 7 r)) ■ S7ip(y)^j dy, (3.31) 



where the Riesz transforms are taken with respect to the y variable. 

The third factor inside the integral is given by Riesz transforms in y of the function 

4>x,r(y) = V'(y) = (u(x + ry) - u(x, r)) ■ Vip(y) = div y ((u(x + ry) - u(x, r))ip(y)\ , (3.32) 

since u is divergence free. The function ip is supported in B\ and we have J B ip(y) dy = 0. Moreover, 
since by assumption I(x,r) < Cf(t)r a for some a £ (0, 1), we have by Theorem 12.41 that \\u(x + ry) — 
u(x, < Cf(t)r a , and therefore || V-'||l° c (Bi) < C f(t)r a , uniformly in x. In fact, ip is also C Q , with 

C a norm bounded by Cf(t)r a (note the scaling in y). 

Gathering these bounds together, we see that RiRj(tp) must be bounded in L°°(M. n ) by Cf(t)r a , uniformly 
in x. Indeed, the Riesz transforms are bounded on C a and on L 2 (since ip is supported on Bi and is bounded 
there, its L 2 norm is also finite), and therefore RiRj{ip) € L 2 {R n ) n C a (R n ) D L°°(R"). 

The Riesz transforms of functions with compact support are not compactly supported. The decay of the 
Riesz transform of a compactly supported function is normally of order — n. However, in this case since the 
function t/j has integral zero, and since it is a derivative, we have that RiRj{xp) decays like |y|~" -2 for \y\ large. 
To see this, let be the Kernel associated with the Riesz transform and, using that div u = 0, we compute 



dz 



R ij{^){v) = J 4 } ( z ) K ij(y - z)dz = J ip(z)(u(x + rz) - u(x,r)) ■ VK i:j (y - z) 

tp(z)(u(x + rz) - u(x,r)) ■ (vKij(y - z) - VKij(y)j dz 

(3.33) 
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in principal value sense, since J ip(u — u) = 0. Letting y be such that \y\ > 2, we obtain that \V 2 Kij{£)\ < 
G\y\ n+2 , for all £ that lies between y and y — z, where C is a sufficiently large dimensional constant, gives the 
desired decay in \y\. 

To summarize, we have proved that 



jCf(t)r a if|y|<2 

S Cf(t)r" 
{ |y|" + 2 



RiRjKujx + ry) - u(x,r)) ■ V<p(y)} < \ c J f ^ if ^ >2 (3-34) 



where C is a universal constant that does not depend on x. Recall that by our assumption and Campanato's 
theorem we also have (u(x + ry) — u(x, r)) < Cf(t)(ry) a . Therefore, we obtain from (13.311 and (I3.34t the 
estimate 



V<C!{t) 2 r 2a - Y [\ \b(x + ry)-b(x,r)\dy+ f \b(x + ry) -b(x,r)\ * dy 

\Jb 2 Jm"\b 2 w\ 



(3.35) 



It is here where it is necessary to make the distinction between the Holder and Morrey-Camapanato cases for the 
a priori assumptions on b, by making the specific choices for b. 

The Holder case 

If b(x, r) = b(x), it is clear that 

\b(x + ry) - b(x) \ dy < CM(x,r/2) < Cg{t)r p (3.36) 

for j3 <S (0, 1]. To bound the tail of the integral arising from the Riesz transforms, we first change variables 
ry = z so that 



i , 2 _ a r i 



\b(x + ry) - b(x)\-—^ dy = r^ / \b(x + z) - b(x)\ 

R"\B 2 \y\ JR"\B 2t 



dz 



\b(x + z)- b(x)\ dcr(z) 1 dp 

2r JdB p P 

- t2 ~° r i (x, |4(i + z) - i(x)] a " {z) ) ip 

<Cr 2 -/ {p n M (x,p))—^dp 

Jlr P 
poo i 

< Cr 2 - a g(t) / -j-^ < Cg(ty. (3.37) 

From d335l . ( T3~36] l, and COT} , we obtain 

V < Cf(t) 2 g(t)r 2a+fi ' 1 . (3.38) 

The BMO case 

In this case we have b(x, r) = (l/r n ) J B ^ b(z) dz. The difference with the Holder case lies in the tail of the 
integral due to the Riesz transform. We use the following classical fact about BMO functions: the difference 
between the mean on B r an B\ r is bounded by 2™ m(l + A) times the BMO norm, for all A > 1. We split the 
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\b{x + ry)-b{x,r)\ 2 _ a dy 



integral in dyadic cylinders to find a bound for each and sum. 

/ \Kx + ry)-b(x,r)\ 1 dy = f 

JS."\B 2 \y\ k=1 J 'B 2k+1 \B 2 

oo 

< E / (M* + r y) - ^ 2fc+v )i + i 5 ^ 2fc+v ) - ^ r )D ok(n+2- a) d y 

fc=l JB 2k + 1 \B 2k 
°° 1 1 /■ 

^E^ylHI^ + ^^y/ |6(x 1 2^r)-6(x,r)|d» 



00 



1 1 

< E ^) IHIbmo + ^T^2 fc « log(l + 2^)\\b\\ BMO < C\\b\\ BM0 . 

fc=i 

Thus, the tail of the integral in the bound of T 5 ( 13.351 ) is bounded by C||6||bmo < Cg(t), as well as the fist term 
in ( 13.351 l. Therefore, in this case (j3 = 0) we also obtain 



V < Cf{t) 2 g {t)r 2a+ P-\ 

The Morrey-Campanato case 

In this case we have b(x, r ) = 0. The same proof as in the Holder case above, but with 6 = 0, shows that 

V < Cf{tfg{t)r 2a+p - 1 (3.39) 
when sup r>0 r~@M(x, r) < 00, and (3 e [—1, 0). □ 

Once we have estimated the three terms on the right side of (13. 16b . the proof of the main theorem is concluded 
as follows. 

Proof of Theorem U . 71 Since the vector drift b is a priori assumed to lie in M 13 , for some f3 e (—1, 1], with 
[b(;t)] M P < g(t), and g g L 2/(1+l3 \ we obtain from dX20l . (TP8t . and ( T3~38l l the estimate 

f'(t)f(t)r 2a < C*f(t) 2 g{t)r 2a +^ - cj(t) 2 r 2a - 2 (3.40) 

which holds for some sufficiently large constant C*, and some sufficiently small positive constant c*. The above 
estimate implies 

< C*g(ty-i c^ 2 < Cg{t)W+V (3.41) 

for some positive constant C = C*(C*, c*, /?). The last inequality was obtained maximizing the expression with 
respect to r. 

On the other hand, if g <G L 2 ^ 1+l3 \ we can choose / to be the solution to the ODE 

fit) = 2C 5 (i) 2/(1+/3) /(t), 

which contradicts the above inequality and makes it impossible for the Holder modulus of continuity to ever be 
invalidated. Note that the above ODE has the explicit solution 

/(t) = e (/o t W /(1+ « d *)/(0), 
which stays bounded for all t. □ 



13 



4 The endpoint case (3 = 



-1 



For the three-dimensional Navier-Stokes equations, i.e. b = u, obtaining regularity of the solutions in the 
endpoint case u € L°°L 3 is highly non-trivial, and this issue was only settled recently by Iskauriaza, Seregin, 
and Sverak in ifTol (see also ifTTI for the case L 00 !! 1 / 2 ). The case (p, q) = (oo, n) on the Ladyzhenskaya-Foias- 
Prodi-Serrin scale is of particular importance as it is the scaling-critical space for the initial data, and it gives the 
borderline space (on the Lebesgue scale) for constructing solutions via the Picard iteration scheme (so-called 
mild solutions). Note that L 3 C A/ -1 . 

It turns out that for the linear system (11.8) -( 11.10h the proof given above fails in the case b G i°°A/ _1 , as 
the constant C blows up as /? — > — 1. The corresponding result for the Navier-Stokes equations (as in iflOl or 
IfTTI ) relies essentially on the nonlinear structure of the equation. In order to obtain a result in this direction for 
the liner equation ( 11.2b , with the method of this article, we need to impose an extra smallness condition. 

Theorem 4.1. Let T > be arbitrary. Assume that b is a divergence-free vector field in L°°([0, T]; M ), and 
let Uq G L 2 H C a for some a G (0, 1). There exists a positive constant e > such that if for all t G (0, T] there 
exists r* (t) > with 

sup sup sup r / \b(x + ry, t)\ dy < e (4.1) 

E[0,T]i6K»0<r<r,(t) J B x 



and 



thenu(-,t) G C a (R™) for all t G (0,Tj. 



T 1 



r*(t) : 



dt < oo, (4.2) 



Remark 4.2. If H^Il^m- 1 < e > me theorem holds trivially. Additionally, note that for any cj) G L J (M' 3 ), or any 
<f> in the closure of in M , we have 

lim r / \4>(x + ry)\ dy = 0, 
Jb 1 

for all x, and hence ( 14. Il l holds for some r* > 0. Therefore, if b is continuous in time with values in L 3 , or 
piecewise continuous with arbitrarily large jumps, or if all jumps are of size smaller than e/2, the conditions 
( 14.1l >- (14.2t are automatically satisfied, with r* (t) being a sufficiently small constant, proving regularity of the 
solution (see also J6] Theorem 3.1] for a similar result in the critical Besov space). Theorem l4. 1 I states that if the 
drift velocity &(•, t) is nicely behaved at sufficiently small scales r*(i), and if these scales do not go to too fast 
at any point t G [0, T], i.e. r st (s) cannot vanish at t with a rate faster than \Jt — s, then the solution is regular 
until T. Note that with respect to time regularity, the conditions of Theorem [4J] above are stronger than merely 
b G LfL 3 ^ but weaker than C t L 3 x . 

Proof of Theoretn\4l\ From (l3~20i i, (l3~28t . (l3~36l l, and ( T3~37l i, we obtain 

2f(t)f(t)r 2a < Cy- 1 f(t) 2 M(x,t,r) - f {t) 2 r 2a ~ 2 

+ C*f(t) 2 r 2a - 1 (M(x,t,r/2) + r 2 - a J™ M(x,t,p)^d p y (4.3) 

where Ai(x, t, r) = J B ^ \b(x + ry, t)\ dy, and by assumption of the theorem we have 

sup rM(x, i, r) < e (4.4) 

0<r<r»(t) 
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and 



suprM(x,t,r) < \\b\\ LaoM -i = B < oo (4.5) 

r>0 * 1 

uniformly in x, t. Inserting the bounds (14 . 41 > — (14 . 5 b into ( 14.31 . we have 

2r 2 -^ log/0) < C*rM(x,t,r) - c, + C*rM(x,t,r/2) 
at 

+ C r r 3 - a (£\ p M(x,t,p))^ dp + J™ (pM(x,t,p))-^ dp*) (4.6) 
for all r > 0. To bound the right side of ( 14.6b . we distinguish the cases r/r* < S, and r/r* > (5, where we let 



l/(3-o) i 



S = ^{[b) '2/- ^ 

Indeed, if r/r* < (5, we have 2?' < r* and B(r /r*) 3 ~ a < BS 3 ~ a < e, so that by ( 14.11 ) implies 

2r 2 ^log/(t) < C*e-c* + 2C*e + C*r 3 - a C a (er a - 3 + Br^ 3 ) 

< 0*(3 + 2C a )e-c*, (4.8) 



for some positive constant C Q > 0, and for all < r < 6 . Therefore, if we choose e as 



0.(3 + 2C a )' 

then we obtain from (14.8b that 



(4.9) 



^log/(t)<0 (4.10) 

for all t € [0, T], and all < r < <5 r* (i). On the other hand, if r > S r* (t), we bound the right side of ( 14.6b by 
making use of (14.5b . namely 

< £ ( 3 C,B + ft Sr - £ ^) < S^^l (4.11, 
for some O a > 0, for all r > <5 r* (i). The proof of the theorem is then concluded since we may choose / as 

which is finite for alii < T by ( |4~2"1 ), and it contradicts ( 14. 10b — ( |4~TTT >. □ 
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